Solving Block Low-Rank Matrix Eigenvalue Problems

نویسندگان

چکیده

To solve large-scale matrix eigenvalue problems (EVPs), a two-step tridiagonalization method using the block Householder transformation (HT) is often employed. Although based on dense arithmetic requires memory storage of O(N2) and an operations O(N3), in this study, these were reduced by approximating low-rank (BLR-) matrices. A special HT for BLR-matrices it are proposed to EVP with real symmetric BLR-matrix. In HT, vectors also formed BLR-matrices. It demonstrated how size m BLR-matrix should be determined confirmed that complexities O(N5/3) O(N7/3), respectively, typical cases when appropriate ∝ N1/3. numerical experiments string free vibration problem known analytical solutions, large eigenvalues, calculated eigenvalues converge toward ones accordance theoretical convergence curves. Owing complexity, was solved about N =300,000, which significantly larger than limit conventional matrices, within reasonable amount time CPU core. For calculation time, faster few tens thousands.

برای دانلود باید عضویت طلایی داشته باشید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Stochastic algorithms for solving structured low-rank matrix approximation problems

In this paper, we investigate the complexity of the numerical construction of the Hankel structured low-rank approximation (HSLRA) problem, and develop a family of algorithms to solve this problem. Briefly, HSLRA is the problem of finding the closest (in some pre-defined norm) rank r approximation of a given Hankel matrix, which is also of Hankel structure. We demonstrate that finding optimal s...

متن کامل

A Density Matrix-based Algorithm for Solving Eigenvalue Problems

A new numerical algorithm for solving the symmetric eigenvalue problem is presented. The technique deviates fundamentally from the traditional Krylov subspace iteration based techniques (Arnoldi and Lanczos algorithms) or other Davidson-Jacobi techniques, and takes its inspiration from the contour integration and density matrix representation in quantum mechanics. It will be shown that this new...

متن کامل

Solving Overdetermined Eigenvalue Problems

We propose a new interpretation of the generalized overdetermined eigenvalue problem (A− λB)v ≈ 0 for two m × n (m > n) matrices A and B, its stability analysis, and an efficient algorithm for solving it. Usually, the matrix pencil {A− λB} does not have any rank deficient member. Therefore we aim to compute λ for which A − λB is as close as possible to rank deficient; i.e., we search for λ that...

متن کامل

A geometric method for eigenvalue problems with low-rank perturbations

We consider the problem of finding the spectrum of an operator taking the form of a low-rank (rank one or two) non-normal perturbation of a well-understood operator, motivated by a number of problems of applied interest which take this form. We use the fact that the system is a low-rank perturbation of a solved problem, together with a simple idea of classical differential geometry (the envelop...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

ژورنال

عنوان ژورنال: Journal of information processing

سال: 2022

ISSN: ['0387-6101']

DOI: https://doi.org/10.2197/ipsjjip.30.538